Drift, Drag and Brownian motion in the Davies-Unruh bath 
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An interesting feature of the Davies-Unruh effect is that an uniformly accelerated observer sees 
an isotropic thermal spectrum of particles even though there is a preferred direction in this context, 
determined by the direction of the acceleration g. We investigate the thermal fluctuations in the 
Unruh bath by studying the Brownian motion of particles in the bath, especially as regards to 
isotropy. We find that the thermal fluctuations are anisotropic and induce different frictional drag 
forces on the Brownian particle depending on whether it has a drift velocity along the direction of 
acceleration g or in a direction transverse to it. Using the fluctuation-dissipation theorem, we argue 
that this anisotropy arises due to quantum correlations in the fluctuations at large correlation time 
scales. 



In 1976, Unruh discovered [H, 0] that a uniformly ac- 
celerated observer, moving in a fiat spacctime, perceives 
the inertial vacuum state to be a thermal bath with a 
temperature T related to the magnitude of its acceler- 
ation g as T = g/2ir. The spectrum of particles seen 
by the observer is known to be isotropic in spite of the 
fact that the uniformly accelerated observer obviously 
has a preferred direction, viz. the direction of acceler- 
ation. Alternatively, one can say that the spectrum of 
quantum fluctuations evaluated along the integral curves 
of the Lorentz boost killing vector are isotropic. Given 
this result, it is natural to ask whether small deviations in 
motion away from the integral curves of the Killing vec- 
tor break this isotropy. Physically, this would correspond 
to the following question: If one moves with, say a uni- 
form velocity with respect to the thermal bath itself, does 
the spectrum of particles as seen by the drifting observer 
depend on the relative orientations of the acceleration 
with respect to the drift velocity? [Motion relative to 
an isotropic spectrum will always introduce the standard 
Doppler anisotropy with a cos 9 dependence with respect 
to the direction of the velocity, even in the absence of 
acceleration etc. This is not what we are interested in. 
We would like to know whether the spectrum depends 
on the relative orientation of velocity and acceleration in 
addition to this standard Doppler shift dependence on 
the direction of velocity alone.] 

We will show in this paper that the spectrum of par- 
ticles seen by the drifting observer does indeed depend 
on whether she is moving along g or in a direction trans- 
verse to it. In particular, the spectrum in the two cases 
match with each other only in the high frequency regime 
while they differ significantly at low frequencies. 

These anisotropics also reflect in the Brownian motion 
of a particle in the Davies-Unruh thermal bath which can 
be seen as follows: The Brownian motion of a particle in 
a thermal bath is described by the Langevin equation 

dv _ 



7(i - t')v{t')dt' + F{t) + F ei 



(1) 



where j(t—t') is the dissipation function, J-(t) is the ran- 
dom part of the force per unit mass due the thermal fluc- 
tuations and F ext is the external force per unit mass act- 
ing on the particle. We calculate the spectrum of parti- 
cles seen by a drifting observer in the Davies-Unruh ther- 
mal bath and show that the observer experiences a drag 
force of the form F = -~jv opposite to its direction of 
drift . This allows us to read-off 7(0;) which is the spectral 
decomposition of the dissipation function j(t—t') appear- 
ing in Eq.([T]) obtained by a Fourier transform. Then the 
result derived below implies that 7g(£ — t') ^ 7 g - J _ (t — t') 
at large correlation time scales. This anisotropy is purely 
quantum mechanical in origin. Consequences of this ef- 
fect on Brownian motion are discussed in the last section. 

We will now proceed to describe the setup we use for 
the calculation. We consider a two level quantum system 
linearly coupled to a quantum field (known as the Unruh- 
Dewitt detector 0]) and analyze its excitation probability 
when it is in motion. It is well known that using linear 
perturbation theory the angular dependence of the ex- 
citation probability can be expressed in the form of an 
integral given by [J] 



dR(E) 
dQ 



drdr 1 e - iE( - T - T ^ 



[(T r -T T ,)-n-(X T -X T , 



■(2) 
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where T T — T(r) and X T = X(r) are evaluated on the 
trajectory of the detector. The angular frequency re- 
sponse, dR(E)/dCl, determines the number density of 
particles emitted by the detector in a given solid angle 
dfl. It has been shown that this direction sensitive detec- 
tor correctly reproduces the isotropic nature of the Unruh 
bath when viewed from the rest frame of the bath [3| and 
hence is ideally suited for our present purpose. We next 
choose a suitable trajectory for the detector such that it 
matches with the mean trajectory of the Brownian parti- 
cle in a thermal bath and calculate the angular response 
for the detector to obtain 'y(E). (Note that since j(t — t') 
is an intrinsic property of the thermal bath, it is inde- 
pendent of the detailed model of the Brownian particle 
chosen.) We will now proceed to calculate *y(t — t') first 
for the case of a detector drifting parallel to the direction 
of acceleration g of the Rindlcr observer. 
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A. Detector with a drift parallel to acceleration 

To study the Brownian motion of a particle in a ther- 
mal bath with a drift along the direction of acceleration, 
we choose a trajectory such that the Rindler observer, 
who is in the rest frame of the bath, perceives this detec- 
tor as moving with a relative velocity v in the direction 
of acceleration, say the x -direction. The trajectory in 
the Rindler frame is then 



vjt' , t = -ft' 



(3) 



where v is a constant and defines the the relative velocity 
of the Brownian particle with respect to the rest frame 
of the Davies-Unruh bath. A suitable coordinate trans- 
formation between the static Rindler observer and the 
drifting observer is x = 7(2;' + vt') and t = 7(4' + vx'). 
In these coordinates, the trajectory is simply x' = and 
the metric becomes 



ds 2 = e 2 9^' e 29ivt' 



[-dt' 2 + dx' 2 } + dy 2 + dz 2 (4) 



From the metric, it is clear that the time co-ordinate 
t' is not the proper time r of the drifting observer but 
is related to the proper time r through r = — (1 — 
e aivt ^(gjv). For our purpose, it is enough to consider 
drift velocities which are very small in magnitude and 
we will work in the regime of linear v assuming u « 1. 
Before proceeding further let us first quantify how small 
a v we need to consider. Let us assume that the detec- 
tor coupling to be switched-on during the time interval 
(— L, L) where L is a very large but finite number such 
that EL ^> 1 holds. This condition basically ensures that 
the contour integration over t — t' in Eq. ([2]) vanishes over 
semi-circle in the lower half of the complex u plane. We 
then choose the velocity parameter v such that gvL <C 1. 
These conditions lead to v <C E/g which tell us that the 
low frequencies (i.e., low compared to g) probed by the 
detector, in this approximation, are limited by choice of 
the parameter v. However, we can still choose a non-zero 
v such that d <C £/g < 1 which allows us to probe the 
effects in the low frequency regime. Under this approxi- 
mation, the proper time r to linear order in v is 



a + g -f 2 



(5) 



Using the above relation, the magnitude of the proper 
acceleration for the trajectory in Eq.© can be found, 
upto linear order in v, to be 



g(r) = \Z~aTa~t = g(l + gvr) 



(6) 



The time dependence of acceleration can be understood 
by noticing that the observer following the trajectory 
given by Eq.© crosses different x = constant surfaces 
at different times. This corresponds to crossing different 
uniformly accelerated hyperbolic trajectories at different 
times, which is essentially described by Eq.®. Further, 



the metric in Eq.© can now be expressed in terms of 
the proper time and we get 

ds 2 = e 2 ^ x ' [-dr 2 + (I + gvT)dx' 2 ] + dy 2 + dz 2 (7) 

(It is interesting to note that metric of the 2-dimcnsional 
x' — t hypersurface is conformal to a 2-dimensional FRW 
metric of an expanding universe with scale factor equal to 
g(r)/g.) We will later show that the factor g(r)/g acts as 
the blue-shifting (or red-shifting depending on the sign 
of v)factor for the temperature as seen by the drifting 
detector. This blueshift in temperature accounts for the 
time dependent particle production which we expect due 
to the non-static nature of the above metric. 

We will now proceed to calculate the angular response 
for the detector using the expression in Eq.®. In Eq.®, 
the integration variables arc in terms of the proper time 
t. However, it turns out to be easier (algebraically) to 
perform the contour integral using the time co-ordinate 
t' in Eq.©. We use Eq.© to make a change of variable 
in Eq.© and obtain: 



dR(E) 
dQ 



ds du (1 + 2gvs) e 



-iEu[l+gvs] 



[T — T' — cos 9 (X — X') - ie] 



(8) 



where we have defined u = (t' x — t'^) and 2s = {t'i+t' 2 ). To 
evaluate the above integral we write the trajectory given 
in Eq.Q in terms of the inertial coordinates expressed 
as function of t' as 



T(0 



(i + gvt) 
9 



sinh gt! ; X 



(i + gyf) 
9 



cosh gt 1 



Substituting the above expressions in Eq.© and expand- 
ing the integrand upto linear order in v results in an inte- 
gral which involves closing the contour in the lower half 
complex u plane and evaluating the residues at the poles 
of the integrand. We assume that the linear order result 
can be obtained by first expanding the integrand to order 
v and then evaluating the residues. The response rate of 
the detector is then given by 



dR(E) 



dn 



(g 2 /4)F 2 (9,s)e- lhu du dR v (E) 



[sinh(f) 



dn 



where 



dR v (E) 
dQ 



(g 2 /4)e~ lEu vucosh(^)du 
F- 3 (9,s)F(9,s) [sinh(^) -ie 



(9) 



(10) 



and 



F(9,s) = [coshgs — cos 9 sinh gs] 
F(8,s) = [cos 9 cosh gs — sinh gs]~ 



(11) 



Here, E = E [1 + gvs) and the overhead (■) denotes dif- 
ferentiation w.r.t s. The first term in Eq.© in the Taylor 
expansion is the familiar term one usually gets while cal- 
culating the response of a detector moving on the Rindler 
trajectory, except that we get E instead of E in the ex- 
ponential. Evaluating the integral in the first term then 
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leads to the Planckian distribution E / (exp f3(s)E — 1) of 
particles with /3 _1 (s) = g(s)/2w. The time dependence 
of the temperature is easily understood in this case as 
being due to the gravitational redshift when the detector 
crosses different x = constant surfaces at different times. 
As is well known, the thermal spectrum of particles seen 
by the Rindlcr observer is isotropic, hence we expect the 
first term also to be isotropic since it differs from the 
former case only in the time dependence of temperature 
/? . This will be evident from our analysis below. 

The second term, which is linear in the velocity, turns 
out to be anisotropic in the spectrum and leads to the 
drag force on the Brownian particle. To calculate the 
required response function, we begin by first performing 
the contour integral over u. Due to the Taylor expansion, 
this turns out to be a simple procedure of evaluating the 
residue of the integrand at the same poles one gets, but 
now of one order higher, while evaluating the first term 
in Eq.©. The response function is then obtained to be 



dR v {E) 3 lx (E) 



dn 



4tt E 



vF 3 (6,s)F- 1 (6,s) (12) 



where, 



7*(E) 



TT 2 f3Eh 



-BE 



1 



2 

Je 



2e 



-BE 



2ir 2 



(1 - e-P E f [' PE ' (3E fi 2 E 2 

(13) 

Hereafter, for notational simplicity, we suppress the time- 
dependence in f3(s). We will show below that jxiE) 
is the spectral decomposition of the dissipation func- 
tion discussed after Eq. ([!]). Note that the angular de- 
pendence of the response function in above Eq. (|12l) is 
in terms of the solid angle dQ defined in the inertial 
frame. One would have expected the angular response 
to be time-independent (except for the time dependence 
in f3(s)) whereas from the F(9, s) dependence, it appears 
as though it depends on the proper time of the trajec- 
tory in a complicated manner. But recall that the natural 
frame of reference to consider when studying Brownian 
motion would be the rest frame of the thermal bath it- 
self which, in the present case, is the Rindlcr frame. As 
we will show, all the peculiar behaviour goes away, when 
we view it in the Rindlcr observer's frame. The relevant 
transformation between the angular co-ordinates of the 
Minkowski frame (9, 4>) and the Rindler frame (9r, </>#) 
is given by 

dil sin6 
dilR = -5-7- -rr^; sin0R 



Jr(1 - vrcosI 



7^(1 - vr cos( 



cos 9 



Luii» - Vr 

cos0r = (14 

(1 — Vr COSV) 

where vr = tanh gs is the relative co-ordinate velocity of 
the Rindler observer with respect to the inertial observer 



and 7h = 1/(1 - v 2 R ). 
F{6,s) = 7^(1 - vr cos 
and hence we get, 



From Eq. tfTTj) . we can note that 
and F(9, s) = -y R (cos9 - v R ) 



dRy(E) = 3 lx (E) 
dVL R 4tt E 



-VCOS Vr 



(15) 



The right hand side of the above equation tells us the 
number of particles of energy E hitting the detector per 
unit time in a given direction vi(9r, 4>r). The force F„ 
from the —n direction on the detector due to this bom- 
bardment is then E(dR v (E) / 'dClfi). So, the net force on 
the Brownian particle will be: 



-1x{E)v 



(16) 



Hence, ^xiE) defined in Eq. ([T3)) is the spectral decom- 
position of the required dissipation function 7 x (i — t'). 

We will next calculate j(t — t') for a detector drifting 
transverse to the direction of acceleration g of the Rindlcr 
observer. 



B. Detector with a drift transverse to the 
acceleration 

To study this case, we choose a trajectory such that 
the Rindler observer will perceive the detector as moving 
with a relative velocity v in the direction transverse to 
the acceleration, say the y - direction. The trajectory in 
the Rindler frame is then 







t 



7r 



y = vjt 



(17) 



where the parameter r is the proper time of the observer 
moving on the above trajectory. 

To calculate the angular response for the detector drift- 
ing along the y axis we proceed in a manner similar to the 
case of the x— drift. However, as we will see, the angular 
response for the y— drift is relatively easier to calculate 
due to the fact that the t' co-ordinate is the proper time 
parameter in this case. We begin by expressing the tra- 
jectory Eq. fTTj) in the inertial coordinates as 

T = — sinh g^T ; X = — cosh g^T ; Y = v^t (18) 
9 9 

Substituting in Eq.([5]) and expanding upto linear order 
in d, we get 



dk{E) 



dCl 



(g 2 /4)F 2 (9,s)e- lIlu du dR v (E) 



[sinh(f ) 



dn 



where 



dR v (E) 
dQ 



vu sin v cos < 



F- 3 (6,s) [sinh( 



2 ' 



(19) 



(20) 



As in the x drift case, the velocity independent term in 
Eq. (|T9|) gives the usual isotropic Planckian distribution 
-©/(exp j3 E — 1) of particles as seen by a detector moving 
on the Rindler trajectory and will not lead to a drag force 
on the detector. The second term linear in the velocity 
is anisotropic and leads to the required drag force. Per- 
forming the contour integral over it, the response function 
then becomes 



dR v {E) _ 3 ly {E) 



dQ 



4tt E 



(vsm9cos(l))F 3 (9,s) (21) 
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whore. 

7 JE) = 5- 

V 3(l + e^) 2 



2 

Je 



2e-? E 3ir 2 1 



(22) 



(The angular response for the detector drifting along the 
y— axis has been obtained before in [B[ in the context of 
an accelerating braneworld model.) We will now proceed 
to calculate the drag force and the dissipation function 
and compare it with the one we obtained in the case of 
the x— drift. We first write the response function in 
terms of the angular co-ordinates of the Rindlcr observer 
using Eq.((T4]) and Eq.(JTTl)- We get: 



dR v (E) = 3n ly {E) 
dfl R ~ 2 E 3 



(usin0RCOS0ii) (23) 



Hence, the net drag force on the detector in the y direc- 
tion is 



F„ 



-ly{E)i 



(24) 



Hence, "j y {E) defined in Eq. (|22"l) is the spectral decom- 
position of the required dissipation function j y (t — t'). 



C. Discussion and Conclusions 

The derivation presented in this letter provides a calcu- 
lation of the dissipation function ^(t — t 1 ) for the Davies- 
Unruh bath from first principles. Surprisingly, one finds 
that 7(t — t') is anisotropic and depends whether one 
is moving in the direction of the acceleration or trans- 
verse to it. The Davies-Unruh effect itself arises due 
to the difference in the pattern of vacuum fluctuations 
seen by inertial and accelerated observers. Therefore we 
can interpret this anisotropy as due to the dependence of 
the fluctuation pattern on both velocity and acceleration 
which will depend on the relative orientation between 
the two. This interpretation is strengthened by the fact 
that Eq. p^)) and Eq. (f2"2")) . for the dissipation function -f x 
and "f y respectively, agree in the high frequency regime 
in which /3E 3> 1. This feature is expected (and a cross- 
check on the calculation) because in the high frequency 
regime the quantum correlations will become more and 
more local and the difference in orientation will have less 
effect. 

In fact, no anisotropy of the type we are discussing will 
arise if we analyze the Brownian drift of a particle in a 
box containing massless particles in thermal equilibrium 
in a constant gravitational field, say, along the x direc- 
tion. In this case the equilibrium temperature of the gas 
will vary along x axis in accordance with the Tolman Q 
relation. An observer, drifting along the direction trans- 
verse to g, will see the spectrum to be Doppler shifted as 
a consequence of the standard Doppler shift in frequen- 
cies Q and will not experience the variation of tempera- 
ture due to Tolman factor. An observer drifting along the 
direction of g will experience, in addition to this standard 



Doppler shift, an extra gravitational redshift in the fre- 
quencies due to her motion along the gravitational field. 
However, the gravitational redshift factor will only make 
the temperature /3 -1 (a;(£)) (obeying the Tolman relation) 
to be time dependent and a function of the position of the 
observer along x which is similar to the time dependence 
we found in the case of parallel drift in the Davies-Unruh 
bath. It is easy to see that the form of the frequency 
response as seen by both the drifting observers would be 
the same and isotropic in the sense that it depends only 
on the direction of v and not on its relative orientation 
with respect to g. (Of course, the spectrum will vary 
with a cos 9 dependence with respect to the direction of 
the velocity in both cases but this is not the anisotropy 
we are interested in.) This is because such an analysis 
is essentially 'classical' in the sense that it can be repro- 
duced by assuming the existence of a bunch of particles in 
a region with a particular classical distribution function 
and doing just the kinematics. No quantum correlations 
come into play. Clearly our result agrees with the above 
analysis and the anisotropy vanishes in this limit. 

As one moves away from the high frequency end of 
the spectrum, the responses of the two detectors start 
to differ significantly. When j3E < 1, we have 7^ 







-2 P -/3E 



/(/3E) while 7 y - P~ 2 (PE)e-^ E . Thus the en- 



ergy dissipated by a Brownian particle due to friction is 
greater by a fraction 1/((3E) 2 because of its mean ve- 
locity in the x direction than the corresponding energy 
dissipated in the y direction. 

Let us consider this result from the perspective of the 
fluctuation-dissipation theorem [f| which relates the dis- 
sipation function to the two point correlation function of 
the random force T at different times as 



7H = ^ PjmHt+t))* 



-itut 



d£ (25) 



In the present context, since we have determined the dis- 
sipation function from first principles, we can now deter- 
mine the characteristics of the fluctuation (or the noise) 
K(u) = (J-(t)F(t + u)) intrinsic to the system using the 
above equation. We find that, at small {{u/ 0) —> 0) cor- 
relation times — which corresponds to high frequency 
range — we get K x rj K y . On the other hand, for very 
long (u) correlation times — which corresponds to the 
low frequency range — both K x and K y vanishes with 
K x {u) / K y {u) oc u 2 because of the different E scaling in 
and 7 U . The behaviour K(u) — > for u — > 00 sig- 
nifies that a Brownian particle which is in equilibrium 
with the thermal bath loses the memory of its path his- 
tory (or its initial conditions) as time progresses. But 
our calculations also show that even though the correla- 
tions do vanish at large correlation times, there are still 
quantum correlations present in the case of x drift which 
go towards zero at a slower rate as compared to the y 
drift. 

We know that the Unruh effect is closely related to 
the Hawking radiation of a non-extremal black hole hori- 
zon, cosmological de-Sitter horizon, etc. Hence, we would 
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expect similar anisotropies to be present in the thermal 
fluctuations in the thermal bath near these horizons too. 
Particularly, the trajectories of particles near the horizon 
would be affected differently for different orientations of 
acceleration and velocity leading to anisotropic dissipa- 
tion effects. It will be interesting to investigate conse- 



quences of this feature. 
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